Splitting of Surface Plasmon Frequencies of Metal Particles in a 

Nematic Liquid Crystal 



Sung Yong Park and D. Stroud 
Department of Physics, The Ohio State University, Columbus, Ohio 43210 

Abstract 



We calculate the effective dielectric function for a suspension of small 
metallic particles immersed in a nematic liquid crystal (NLC) host. For a 
random suspension of such particles in the dilute limit, we calculate the effec- 
tive dielectric tensor exactly and show that the surface plasmon (SP)resonance 
of such particles splits into two resonances, polarized parallel and perpendicu- 
lar to the NLC director. At higher concentrations, we calculate this splitting 
using a generalized Maxwell-Garnett approximation, which can also be ap- 
plied to a small metal particle coated with NLC. To confirm the accuracy of 
the MGA for NLC-coated spheres, we also use the Discrete Dipole Approx- 
imation. The calculated splitting is comparable to that observed in recent 
experiments on NLC-coated small metal particles. 
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An important goal in optics is to prepare materials with optical properties easily con- 
trollable by external parameters such as temperature or electromagnetic fields. Nematic 
liquid crystals (NLC's) are useful materials of this kind [1-4]: the optical axis of an NLC, 
known as the director, can be rotated by a weak applied electric field thereby altering 
its transmission and refiection coefficients. Several authors have suggested that such control 
might be achieved by incorporating NLC's into composite media which are structured on a 
scale comparable to the wavelength [5-11]. 

Recently, there have been new optical applications of nanoparticles on a scale smaller 
than the wavelength [12]. Among the applications are nanolenses [13], plasmonic nanoan- 
tennas [14], and fight energy transfer through linear chains of nanoparticles [15-17]. Thus, 
it would be of great interest if the use of NLC's to control fight propagation could be ex- 
tended to systems with structure on this scale. Indeed, a recent experiment has shown that 
the optical properties of NLC-coated metal nanoparticles can be controlled by a dc electric 
field [18]. 

In this Letter, as a first step towards a theoretical treatment of such systems, we present 
an exact calculation of the optical properties of a dilute suspension of metallic nanoparticles 
in an NLC host. We find that the usual surface plasmon (SP) resonance of the metal 
nanoparticles splits into two resonances, polarized parallel and perpendicular to the NLC 
director. In addition, we develop a Maxwell-Garnett approximation (MCA) which is suitable 
either for higher concentrations of particles in an NLC host or for a single nanoparticle 
coated with an NLC. Wc also use the discrete dipole approximation (DDA) to confirm the 
accuracy of the MGA for NLC-coated spheres. For typical NLC parameters, we find that the 
calculated SP splitting for coated spheres is large enough to be experimentally observable in 
both absorption and transmission, and, indeed, is of a magnitude comparable to that found 
in experiments [18]. 

First, we consider the optical properties of a dilute suspension of spherical metallic 
nanoparticles in an NLC host, as illustrated schematically in Fig. 1 (a). The metal particles 
are assumed to have a complex but isotropic dielectric function e^(a;), which we take simply 
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to be of the Drude form, 

ern{uj)^l-ujl/[oj{uj + i/T)l (1) 

where ujp is the plasma frequency and r is a characteristic relaxation time. The NLC host 
is taken to be a uniaxial dielectric with dielectric tensor If^Lc having principal components 
e_L = e_L, and ey = Up assumed real and frequency-independent. If the metal particle 
radius is small compared to either its skin depth or the wavelength A, the composite optical 
properties are well described by a complex effective dielectric tensor, leioj). We also assume 
that the NLC is homogeneous, with principal axes which point in the same directions at 
every point in the space it occupies. (The limitations of this assumption are discussed briefly 
below.) With these approximations, 65(0;) is also a tensor with principal components ee,_L(a;), 
ee,±(a;), and ee,||(a;). 

We calculate ee{oj) in the limit that the volume fraction of inclusions p <^1. In actuality, 
ee for this geometry has been calculated in Ref. [19] for the formally analogous case of 
conducting particle in a host with an anisotropic d. c. conductivity. The result for the 
dielectric function is 

eeH = ei+p[i-(5ef]-^5e, (2) 

where 5e — €2 — ii is the difference between the diagonal dielectric tensors 62 and ei of the 
inclusion and the host, and F is a 3 x 3 "depolarization tensor" which is diagonal in the 
same frame of reference as i^. In the case of Fig. 1(a), ei corresponds to Inlc and 62 to e^l, 
where 1 is a 3 x 3 unit matrix. Equation (2) is a matrix equation for q but because both 
F and 6e are simultaneously diagonalizable, it is also diagonal and reduces to three separate 
algebraic equations for the three components of eg. If the metallic inclusions are spherical 
particles, then F has principal components Fj^, Fj^, and Fy, where 

F|| = -(1 - yr^sin-^ V9/V9)/{Hg), (3) 

and 
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r± = -^(r,! + sin-i ^l^g^el), (4) 

where g —\ — ex/e||. 

The positions a;a (a = || or ±) of the two peaks can be obtained analytically from the 
equation 

1 - VJe^ = 0, (5) 

where a = || or ±. For a Drude dielectric function in the limit lo^t — > oo, this equation 
reduces to 

^oV^p = -ra/[l-(l-6a)r„], (6) 

where is given by Eqs. (3) and (4). To present illustrative results, we assume that the host 
is the NLC known as E7, which has principal indices of refraction = = 1.52, n\\ = 
= 1.75, as used in the experiments in Ref. [18]. For these n± and n\\, ujQ^a/^p — 0.4068 
and 0.4022. Here the difference is a small but observable amount, and is robust against small 
changes of either index of refraction. It is of interest to compare the calculated splitting 
with that observed in the experiments. For p = 0.01, our calculated sphtting between the || 
and _L SP frequencies is about 1%. By comparison, that observed in Ref. [18] appears to be 
around 3.4%. 

At higher concentrations, can only be computed approximately, e. g., using the 
Maxwell-Garnett approximation (MGA). The MGA is beheved to be best suited to a geom- 
etry in which the inclusions are preferentially surrounded by the host (NLC, in this case). 
If the director fi of the NLC is independent of position, then the MGA takes the form 



^MGA = 62 + P i - (1 - p)V5e 
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51. 

Next, we consider the optical properties of a dilute suspension of coated particles. The 
coating is assumed to consist of a thin layer of NLC. The geometry in this case is shown 
schematically in Fig. 1(b). We may calculate the effective dielectric function of this suspen- 
sion in two stages. First, we calculate the dielectric function e^oat of the coated particle using 
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the MGA, and assuming that the director n in the coating is uniformly oriented. Thus, we 
take 

^coat = ^NLC + p' I - {I - p')rSe ^ Se, (7) 

where s is the thickness of the coating, a is the radius of the included metallic particle, 
p' — a? / (a + s)^ is the volume fraction of metal in the coated particle, and 6e — e^l — e^vLC- 
In the second stage, we assume that a volume fraction p <C 1 of coated particles is embedded 
in an air host, and calculate the effective dielectric tensor Ce of this suspension using eq. (2). 
To carry out this latter calculation, we take ei — l,5e — icoat — 1, and use F = —1/3, the 
depolarization tensor for a spherical particle in an isotropic host of unit dielectric constant. 
The resulting q is diagonal with nonzero components 

^e,i, = 1 + Pt^, (8) 

where 

f _ ^coat,n — 1 /qv 

l + (l/3)(ee.a*,^-l)' 

and the subscript ix — || or _L denotes a component in the direction parallel or perpendic- 
ular to h. By analogy with Eq. (5), we can determine the positions uj^ of the two peaks 
analytically from the equation 

l + (l/3)(eeo„t,^-l)=0, (10) 

where /i = || or ±. 

We have used Eq. (10) to obtain the SP resonance frequencies of a dilute suspension of 
NLC-coated metal particles. We consider two different angles between the director and the 
polarization of incident light, namely 0" and 90°. We use the Drude dielectric function of Eq. 
(1) for e^(a;), in the limit ujpT oo, and the dielectric tensor of E7 for that of the coating. 
The resulting frequencies of the SP resonances are shown in Fig. 2 as full and dashed hues. 

To see how well this MGA approach works for a dilute suspension of coated metal 
particles, we have compared the results to DDA calculations of the normalized peak position 



of scattered amplitude [20]. In the DDA, the scattering object is represented as a collection 
of point dipoles on a simple cubic grid of points, with n as large as 120; we increase 
n until the peak position of scattered amplitude converges. The magnitudes of the point 
dipoles are chosen to reproduce the scattering properties of the original medium [21]. We use 
the Drude dielectric function of Eq. (1) for em{uj), with cUpT — 7.7, and the dielectric tensor 
of E7 for that of the coating. The peak positions, calculated in the DDA and interpreted as 
SP frequencies, are shown in Fig. 2 as open squares and diamonds. These calculations show 
that, for p' > 0.5, the MGA and DDA results are very close. For a — 30nm this result means 
that for s ^ 8nm, the MGA result is reasonable. Thus, the Maxwell- Garnett approach is 
fairly accurate for the optical response of nanoparticles with a thin NLC coating, when the 
director is uniform in the coating. 

Finally, we discuss the limitation of our present calculations. In a plausible NLC layer 
thickness range, the calculated splitting in Fig. 2 between the SP frequencies for light po- 
larized parallel and perpendicular to n is 1~2%, roughly comparable to experiment [18]. 
However, to make a detailed comparison to experiment, we would need to relax some of the 
assumptions made in the present calculations; the assumption that the NLC is homogeneous, 
near the surface of a colloidal particle, may be very difficult to achieve experimentally, since 
in real experimental situations, the NLC director is substantially perturbed near the surface 
of a colloidal particle. The characteristics of the surface of the particle, which are referred 
to as anchoring conditions, can yield various nontrivial configurations of the nematic liquid 
crystal [22]; these anchoring conditions will vary with the experimental preparation [23]. 
Thus, a more quantitatively accurate theory of the optical properties of suspensions should 
include these effects, though the present approach should provide a qualitatively reason- 
able first step. Calculations which take into account these complications will be presented 
elsewhere [24]. 

This work has been supported by NSF grant DMR 01-04987. We would like to thank 
Dr. S. Maier for useful conversations. Numerical calculations in this paper were carried out 
using the facilities of the Ohio Supercomputer Center. 
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FIGURE CAPTIONS 



FIGURE 1. Schematic diagram of the geometries we consider here, (a) metal nanopar- 
ticle embedded in an NLC host, (b) NLC coated metal nanoparticle in air host. 

FIGURE 2. Calculated surface plasmon frequencies a;o,||(p')/wo and a;o,±(p')/<^o, cor- 
responding to incident light polarized parallel and perpendicular to the NLC director, for 
metal nanoparticles coated with NLC and embedded in an air host, plotted as a function of 
volume fraction p' of the coated particle which is composed of metal. The coating is assumed 
to be the NLC E7, which has indices of refraction = 1.75, n± — 1.52. The metal particles 
are taken to have the Drude dielectric function, and a;o = a;o(p' = 1). Full and dashed lines: 
MGA; squares and circles: DDA. 



9 



FIGURES 




(b) 

FIG. 1. 
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